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PREFACE 

This  report  was  prepared  by  Dr.  Erik  H.  Vanmarcke,  Professor  of  Civil 
Engineering  at  Massachusetts  Institute  of  Technology  (MIT) ,  as  part  of  on¬ 
going  work  at  the  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  in 
Civil  Works  Investigation  Studies,  "Methodologies  for  Selecting  Design  Earth¬ 
quakes,"  sponsored  by  the  Office,  Chief  of  Engineers,  U.  S.  Army.  This  is  the 
fourteenth  of  a  series  of  state-of-the-art  reports  on  methodologies  for  de¬ 
termining  the  severity  of  bedrock  motion  during  earthquakes. 

Preparation  of  the  report  was  under  the  direction  of  Dr.  E.  L. 

Krinitzsky,  Engineering  and  Rock  Mechanics  Division  (EG&RMD) ,  Geotechnical 
Laboratory  (GL) .  General  direction  was  by  Mr.  J.  P.  Sale,  Chief,  GL,  and 
Dr.  D.  C.  Banks,  Chief,  EG&RMD.  Tha  author  wishes  to  express  his  apprecia¬ 
tion  to  Professor  Robert  V.  Whitman  of  MIT  for  his  helpful  advice  during  the 
conduct  of  the  study  and  the  preparation  of  this  report.  Mr.  Paul  Lai  of  MIT 
assisted  in  carrying  out  many  of  the  computations. 

COL  J.  L.  Cannon,  CE,  and  COL  Nelson  P.  Conover,  CE,  were  Directors  of  WE 
during  the  period  of  this  study.  Mr.  F.  R.  Brown  was  Technical  Director. 
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STATE-OF-THE-ART  FOR  ASSESSING  EARTHQUAKE  HAZARDS  IN  THE  UNITED  STATES 

REPRESENTATION  OF  EARTHQUAKE  GROUND  MOTION: 

SCALED  ACCELEROGRAMS  AND  EQUIVALENT  RESPONSE  SPECTRA 

PART  I:  INTRODUCTION 

1»  A  most  important  step  in  seismic  safety  evaluation  of  civil  engineer¬ 

ing  facilities  is  the  specification  of  the  seismic  input.  Depending  on 
the  type  of  structure,  its  importance,  and  the  type  of  analysis  contem¬ 
plated,  different  characterizations  of  the  input  earthquake  motion  are 
appropriate.  In  present  earthquake  engineering  practice,  the  seismic  in¬ 
put  will  likely  be  expressed  either  in  terms  of  a  set  of  smooth  response 
spectra  or  by  a  selection  of  one  or  more  representative  "time  histories" 
of  earthquake  motion. 

2.  At  a  given  site,  different  representations  of  essentially  the  same 
seismic  input  may  be  required  to  investigate  different  classes  of  struc¬ 
tures  or  to  perform  different  types  of  analysis.  For  example,  at  a  dam 
site,  the  engineer  may  want  to  assess  the  seismic  safety  of  an  earth  em¬ 
bankment  by  time  history  analysis,  and  the  concrete  spillway  and  miscel¬ 
laneous  appurtenant  structures  by  the  response  spectrum  method  of  analysis. 
In  a  study  of  seismic  behavior  of  the  earth  embankment  itself  for  stability, 
deformation  and  liquefaction  potential,  the  engineer  might  alternately 
represent  the  earthquake  by  a  pseudo-static  force,  a  response  spectrum  and 

a  time  history. 

3.  This  report  is  concerned  with  the  question  of  compatibility  between 
different  representations  for  essentially  the  same  seismic  input.  But  since 
different  modes  of  specifying  the  seismic  input  are  intimately  related 

to  different  procedures  of  analysis,  it  seems  logical  that  criteria  for 
compatibility  be  applied  to  the  end  product  of  the  seismic  analysis,  i.e., 
to  the  predicted  performance  or  the  predicted  seismic  load  effect.  The  aim 
(of  the  search  for  compatibility)  is  to  permit  achievement  of  consistent 
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levels  of  safety,  Independent  of  the  representation  of  seismic  input  and  of 
the  mode  of  seismic  analysis. 

4.  The  report  examines  the  alternatives  for  specifying  earthquake  ground 
motions  which  serve  as  input  for  seismic  analysis  and  design.  The  focus 
is  not  on  how  the  characteristics  of  the  design  earthquake  (say,  magnitude 
and  source-to-site  distance)  are  selected,  or  on  how  the  maximum  ground 
acceleration  is  obtained.  Rather,  the  report  deals  with  the  characteriza¬ 
tion  of  ground  motions  conditional  on  the  "design"  maximum  ground  accelera¬ 
tion  (or  velocity),  or  conditional  on  magnitude,  distance  and  local  site 
geological  condition. 

5.  The  report  presents  the  results  of  a  major  study  aimed  at  evaluating 
the  consequences  of  the  widespread  practice  of  scaling  accelerograms  and 
smooth  response  spectra  on  maximum  acceleration.  In  the  WES  report 
series  on  "Methodologies  for  Selecting  Design  Earthquakes,"  Krinitzsky 

and  Chang  (1)  propose  that  maximum  motion  levels  (acceleration  or  velocity) 
be  used  as  a  basis  for  rescaling  selected  existing  records  or  for  the  gene¬ 
ration  of  appropriate  synthetic  records.  With  respect  to  limitations  on 
scaling,  they  state:  "Because  of  uncertainties  concerning  the  appropriate¬ 
ness  of  rescaling  any  single  earthquake  record,  several  records  should  be 
used.  However,  strong-motion  records  should  be  selected  that  require  as 
little  rescaling  as  possible.  If  a  record  has  to  be  rescaled  by  as  much 
as  a  factor  of  4,  the  records  should  be  discarded."  The  study  of  the  ef¬ 
fect  of  scale  factors  consists  of  a  series  of  statistical  analyses  of  the 
response  of  simple  systems  to  a  set  of  140  horizontal  components  of  70 
earthquake  records.  Its  purpose  is  to  bring  about  a  better  understanding  of 
the  effects  of  scaling  earthquake  records  or  response  spectra,  and  to  assess 
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the  degree  of  conservatism  inherent  in  design  procedures  that  are  based  on 
the  use  of  scaled  response  spectra  or  scaled  earthquake  records . 

6.  In  parallel  with  the  study  on  the  effect  of  scaling,  and  in  an  attempt 
to  explain  its  results,  the  frequency  content  and  the  strong-motion  dura¬ 
tion  of  the  set  of  140  accelerograms  are  also  examined.  The  dependence 

of  duration  and  dominant  frequency  on  magnitude  and  distance,  and  their 
(indirect)  correlation  with  maximum  acceleration,  are  discussed  in  Part  IV 
of  this  report. 

7.  Significant  biases  are  indeed  found  to  be  in  evidence  in  the  study 

of  the  effect  of  accelerogram  scaling.  These  biases  can  be  avoided  if  the 
earthquake  engineer  uses  smooth  response  spectra  that  are  site-specific, 
i.e.,  compatible  with  the  site  seismic  condition  under  study.  In  pursuit 
of  this  idea,  the  last  part  of  the  report  develops  and  illustrates  a  pro¬ 
cedure  for  obtaining  site-specific  smooth  response  spectra,  it  is  intended 

to  serve  the  needs  of  engineers  who  must  assess  the  seismic  safety  of 
projects  (such  as  dams  and  appurtenant  structures)  for  which  seismic  in¬ 
put  is  required  in  the  form  of  smooth  response  spectra  as  well  as  time  his¬ 
tories  . 

8.  It  is  assumed  that  a  catalog  of  strong-motion  earthquake  records  is 
available,  from  which  the  engineer  selects  a  number  of  time  histories  that 
are,  as  much  as  possible,  representative  of  the  site  seismic  condition  being 
investigated.  The  proposed  procedure  generates  a  set  of  smooth  response 
spectra  that  can  represent  the  site  seismic  condition  as  well  as  the  time 
histories.  The  proposed  procedure  is  illustrated  for  two  site  conditions 
characterized  in  terms  of  magnitude,  source-to-site  distance  and  local 
geology  (soil  vs  rock).  In  the  illustration,  it  is  assumed  that  the  family 
of  available  accelerograms  is  the  set  of  140  accelerograms  mentioned  earlier. 
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PART  II:  REPRESENTATION  OF  EARTHQUAKE  GROUND  MOTION 


9.  From  the  perspective  of  earthquake  engineering,  the  most  important 
parameter  characterizing  ground  motion  is  the  maximum  ground  acceleration. 
Seismic  design  criteria  are  most  often  expressed  in  terms  of  maximum  ac¬ 
celeration.  A  design  value  of  maximum  acceleration  may  be  obtained  from 
an  attenuation  relationship  in  function  of  a  specified  magnitude  and  dis¬ 
tance,  from  correlations  with  a  specified  M.  Mercalli  Intensity,  or  direct¬ 
ly  from  a  regional  seismic  risk  map  which  corresponds  to  an  appropriate 
mean  return  period. 

10.  The  maximum  acceleration  serves  as  an  anchor  for  more  detailed  ground 
motion  descriptions.  The  precise  seismic  input  format  obviously  depends 
on  the  type  of  analysis  to  be  performed.  At  present,  whenever  some  form 
of  dynamic  (as  opposed  to  pseudo-static)  analysis  is  contemplated,  the 
seismic  input  takes  the  form  of  either  time  histories  or  smooth  response 
spectra.  Both  can  be  conveniently  scaled  w.r.t.  maximum  acceleration. 

11.  Time  history  analysis  is  commonly  used  to  evaluate  the  response  of 
earth  structures  (to  check  stability,  deformation,  liquefaction  potential) , 
secondary  systems  (equipment)  attached  to  major  structures,  or  systems  whose 
complex,  often  nonlinear,  behavior  is  under  study. 

An  advantage  of  time  history  analysis  is  that  it  is  quite  generally 
applicable.  Repeated  time  history  calculations  of  response  to  different 
accelerograms  (all  scaled  to  the  same  maximum  acceleration)  can  provide 
important  information  about  variability  of  seismic  response,  but  such  ex¬ 
tensive  calculations  are  often  prohibitively  expensive  in  practice. 

12.  The  response  spectrum  method  of  analysis  yields  response  predic¬ 
tions  for  multidegree  linear  (or  linearized)  systems  by  combining  modal 
responses  that  are  taken  directly  from  a  smooth  response  spectrum.  For 
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an  important  class  of  multidegree  inelastic  systems,  the  response  spec¬ 
trum  method  also  provides  useful  predictions  if  the  (elastic)  response 
spectrum  is  first  modified  by  an  "inelastic  response  ratio"  which  depends 
on  the  natural  period  and  on  the  amount  of  inelastic  action  permitted  in 
the  structure  (i.e.,  on  the  "design"  ductility  factor). 

A  common  procedure  for  generating  design  response  spectra  is  to  mul¬ 
tiply  the  design  maximum  ground  acceleration  by  a  set  of  scaled  (unit  ac¬ 
celeration)  response  spectra.  It  is  well  known  that  the  shape  of  the  scaled 
response  spectra  presently  used  in  seismic  design  of  nuclear  power  plants 
(prescribed  by  NRC  Reg.  Guide  1.60)  was  obtained  from  the  results  of  sta¬ 
tistical  processing  of  the  response  spectra  of  a  large  number  of  accelero¬ 
grams,  all  scaled  to  a  common  maximum  acceleration  (2,3).  The  background 
studies  by  Newmark  and  Blume  (see  Ref.  2)  involved  calculations  of  response 
spectra  of  33  and  28  components  of  horizontal  motion,  respectively.  The 
mean  and  the  mean  plus  one  standard  deviation  of  the  spectral  amplifica¬ 
tion  factors  were  computed.  The  latter  level  was  selected  as  the  level 
appropriate  for  use  in  design. 

A  third  method  of  analysis  which  has  not  yet  found  its  way  into  prac¬ 
tice,  but  which  holds  considerable  promise,  is  based  on  random  vibration. 
This  method  shares  the  qualities  of  simplicity  and  economy  with  the  response 
spectrum  method,  and  it  does  permit  assessment  of  the  variability  of  seis¬ 
mic  response  predictions  (and  the  probability  that  extreme,  or  permissible, 
response  levels  will  be  exceeded) .  The  input  for  a  seismic  random  vibra¬ 
tion  analysis  takes  the  form  of  the  spectral  density  function  and  the  dura¬ 
tion  of  strong  ground  motion.  A  basic  property  of  the  spectral  density 
function  is  that  its  integral  over  all  frequencies  is  equal  to  the  rms 
ground  acceleration.  This  quantity  is  in  turn  proportional  to  the  maximum 


ground  acceleration.  The  point  is  that  the  seismic  input  required  for  ran¬ 
dom  vibration  analysis  can  also  be  scaled  on  peak  acceleration  by  adjusting 

the  area  under  the  spectral  density  function  (4,5). 

15.  Spectral  density  functions  can  serve  as  a  starting  point  for  generat¬ 

ing  artificial  earthquake  time  histories.  Moreover,  all  three  modes  of 
ground  motion  representation  (response  spectra,  spectral  density  functions, 
and  time  histories)  come  into  play  when  time  histories  are  synthesized  based 
on  the  criterion  that  their  corresponding  response  spectra  must  closely 
match  a  set  of  prescribed  smooth  response  spectra  (5) . 
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PART  III:  EFFECTS  OF  ACCELEROGRAM  SCALING 
ON  SEISMIC  PERFORMANCE  PARAMETERS 


Purpose  of  the  Study 

16.  The  purpose  of  the  study  described  in  this  section  of  the  report  is 

to  quantify  errors  that  may  result  from  scaling  of  earthquake  records  or 
response  spectra  on  maximum  acceleration  -  a  procedure  which  is  so  much 
a  part  of  current  earthquake  engineering  practice.  The  idea  is  that  max¬ 
imum  acceleration  is  only  one  of  a  number  of  significant  parameters  char¬ 
acterizing  earthquake  ground  motion.  Other  parameters  that  may  affect 
the  amplitude  and  the  shape  of  response  spectra  (and  more  generally,  the 
seismic  performance  of  structures)  are  the  strong-motion  duration  and  the 
predominant  period  of  the  ground  motion.  To  the  extent  that  the  maximum 
acceleration  is  correlated,  either  directly  or  indirectly,  with  one  of 
these  parameters,  biases  will  be  introduced  by  scaling  response  spectra 
or  accelerograms. 

The  Set  of  Accelerograms  Used  in  the  Analysis 

1  7.  The  studies  were  carried  out  using  a  data  set  of  140  horizontal  com¬ 

ponents  of  70  western  United  States  strong-motion  records  corresponding 
to  different  event/site  pairs.  It  is  the  same  data  set  selected  by  McGuire 
(6,7)  and  McGuire  and  Barnard  (8)  to  avoid  bias  due  to  unusually  large  num¬ 
bers  of  records  from  a  single  event,  i.e.,  the  1971  San  Fernando  earth¬ 
quake.  Pertinent  characteristics  of  these  records  are  listed  in  Table  1. 

18.  A  scattergram  of  magnitudes  and  epicentral  distances  is  shown  in 

Fig.  1.  The  near-field/far-field  demarcation  line,  based  on  the  tabula¬ 
tion  presented  in  the  report  by  Krinitzsky  and  Chang  (9)  is  also  shown. 
Eleven  sites  (22  records)  were  classified  by  McGuire  as  "rock"  sites,  and 
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Table  1  (Continued) 
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59  sites  (118  records)  as  "soil"  sites. 

Seismic  Response  Quantities  Investigated 

19.  The  effect  of  scaling  was  evaluated  in  terms  of  the  response  of  var¬ 
ious  simple  elastic  and  elasto-plastic  systems,  and  also  in  terms  of  the 
number  of  equivalent  cycles  to  liquefaction  failure  in  a  soil  mass.  All 
these  seismic  response  quantities  were  computed  for  each  of  the  140  ac¬ 
celerograms.  Before  describing  specific  statistical  analysis  procedures 
as  well  as  results,  a  brief  review  is  presented  of  the  definitions,  ter¬ 
minology,  and  notation  used  in  referring  to  the  various  response  quantities. 

Response  Spectra 

20.  The  maximum  response  of  a  simple  linear  one-degree  structure  with  a 

natural  frequency  and  damping  ratio  8  may  be  represented  in  terms  of 

relative  displacement  (S„) ,  the  pseudo-velocity  (S„  =  oi  S„) ,  or  pseudo- 

1)  V  n  D 

acceleration  (S^  =  oo^2  S^) .  The  pseudo-velocity  is  approximately  equal 
to  the  maximum  relative  velocity  of  the  one-degree  structure,  while  the 
pseudo-acceleration  is  close  to  its  maximum  absolute  acceleration.  Re¬ 
sponse  spectra  are  plots  of  these  maximum  response  parameters  as  func¬ 
tions  of  natural  frequency,  or  natural  period  (T^  =  2-ir/io^)  ,  for  different 
values  of  damping.  The  response  spectra  computed  for  this  study  are  all 
pseudo-acceleration  spectra. 

Inelastic  Response  Snectra 

21.  Newmark  (10,  11)  developed  a  simplified  procedure  for  determining  the 
relationship  between  the  maximum  response  of  a  single-degree-of -freedom 
elasto-plastic  system  and  the  maximum  response  of  the  associated  elastic 
system.  The  result  is  the  inelastic  response  spectrum  which  consists  of 
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an  "acceleration  spectrum"  and  a  "(relative)  displacement  spectrum."  These 
may  be  plotted  as  a  function  of  the  (initial  undamped)  natural  period  for 
a  specified  level  of  the  ductility  ratio  y  and  the  damping  factor  (3.  The 
ductility  ratio  y  is  defined  as  the  ratio  of  the  maximum  deflection  to 
the  limit  elastic  deflection,  y  =  y  /y  ,  as  shown  in  Fig.  2.  The  inelas- 
tic  acceleration  response  spectrum  will  be  denoted  by  S^(y) . 

22.  Starting  from  the  pseudo-velocity  elastic  response  spectrum  plotted 
on  log-log  paper,  the  inelastic  response  spectrum  is  developed  as  follows. 

In  the  period  range  where  displacement  or  velocity  are  amplified,  the  in¬ 
elastic  displacement  spectrum  is  identical  to  the  elastic  spectrum.  The 
inelastic  acceleration  spectrum  in  these  regions  is  obtained  by  dividing 
the  elastic  spectrum  by  the  ductility  ratio  y.  In  the  amplified  acceler¬ 
ation  range  of  the  response  spectrum,  the  inelastic  acceleration  spec¬ 
trum  is  obtained  by  dividing  the  elastic  spectrum  by  /2y-l  .  At  very  high 
frequencies,  the  elastic  and  inelastic  response  spectra  are  identical. 

The  spectrum  is  completed  by  drawing  a  straight  line  from  the  zero-period 
acceleration  line  to  the  amplified  acceleration  line.  Note  that  the  in¬ 
elastic  displacement  and  acceleration  spectra  always  differ  by  the  factor 
y.  The  ratio  between  the  elastic  and  the  inelastic  acceleration  response 
spectra,  S^/S^(y),  is  referred  to  as  the  inelastic  acceleration  response 
ratio  (see  Fig.  2b) . 

23.  For  example.  Fig.  3  shows  a  smooth  elastic  response  spectrum  and  the 
corresponding  response  spectrum  for  a  damping  of  2  percent  and  a  ductility 
ratio  of  4. 

24.  Requirements  to  resist  damage  for  structures  and  components  are  often 
expressed  in  terms  of  an  allowable  ductility  factor.  The  elastic  condi¬ 
tion  corresponds  to  y  =  1,  while  values  y  >_  5  may  be  appropriate  when  a 
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Figure  3.  Relationship  between  smooth  elastic  and  inelastic  response 
spectra.  The  inelastic  response  spectrum  is  for  a  ductility  factor 


y  = 


considerable  amount  of  energy  can  be  absorbed  in  inelastic  deformation. 
Ductility  factors  appropriate  for  use  in  design  are  discussed  in  Refer¬ 
ences  10  and  11. 

Equivalent  Number  of  Cycles 

25.  Dynamic  tests  for  soil  liquefaction  potential  measure  the  number  of 
cycles  at  constant  shear  stress  amplitude  necessary  for  excessive  strain 
to  develop  in  a  soil  sample.  The  strength  of  the  soil  under  conditions 
of  cyclic  loading  may  be  represented  by  a  plot  of  shear  stress  amplitude 
versus  number  of  cycles  to  cause  initial  liquefaction  (12,  13). 

26.  Shear  stress  time  histories  generated  by  real  earthquakes  consist 
of  a  succession  of  peaks  and  valleys  with  varying  amplitudes.  Lee  and 
Chan  (14)  developed  a  simple  procedure,  which  parallels  the  probabilistic 
Miner  approach  (15)  to  fatigue  damage  of  metals,  to  evaluate  the  equivalent 
number  of  cycles  in  an  acceleration  time  history.  The  result  will  depend 
on  the  stress  level,  or  on  the  specified  fraction  of  the  maximum  accelera¬ 
tion,  for  which  equivalence  is  sought,  and  also  on  the  details  of  the 
"liquefaction  law"  of  the  soil.  The  particular  relationship  used  in  this 
study  expresses  the  ratio  of  "equivalent  acceleration"  to  peak  acceleration 
versus  number  of  cycles  to  cause  initial  liquefaction. 

Statistical  Analysis  Procedure 

27.  The  specific  response  quantities  selected  for  the  statistical  study 
of  the  effect  of  scaling  are  all  defined  in  such  a  way  that  they  have  the 
property  of  being  "invariant  with  respect  to  scaling  of  the  accelerogram." 
This  means  that  scaling  a  particular  accelerogram  (i.e.,  multiplying  all 
its  accelerations  by  a  common  scale  factor)  will  not  change  the  calculated 
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values  (for  that  accelerogram) .  The  following  response  quantities  have 
this  property  of  invariance: 

a)  Acceleration  amplification  factor  =  ratio  of  the  acceleration 
response  spectrum  to  the  maximum  acceleration.  This  factor  can 
be  evaluated  for  either  elastic  or  inelastic  response  spectra. 

It  depends  on  natural  period  T^,  damping  ratio  3,  and  ductility 
factor  p.  The  value  p=  1  implies  elastic  response  spectra. 

b)  Inelastic  acceleration  response  ratio  =  ratio  of  the  ordinate 
of  the  elastic  response  spectrum  to  that  of  the  inelastic  response 
spectrum  (for  the  same  natural  period  and  damping,  and  for  a  given 
value  of  the  ductility  factor) . 

c)  Equivalent  number  of  cycles  corresponding  to  an  acceleration 
level  that  is  a  specified  fraction  of  the  peak  acceleration.  The 
invariance  property  rests  on  the  assumption  that  the  "liquefaction 
laws"  can  be  represented  by  a  plot  of  a/ a^  versus  N  (in  which  a/ap  =  1 
corresponds  to  N  =  1,  and  N  increases  as  a/ap  decreases). 

28.  Given  a  set  of  accelerograms,  it  is  possible  to  calculate  the  average, 

the  standard  deviation,  and  to  construct  a  histogram,  for  any  of  the  re¬ 
sponse  measures  just  mentioned.  [Precisely  this  kind  of  statistical  analy¬ 
sis  (for  the  acceleration  amplification  factor)  gave  rise  to  the  develop¬ 
ment  of  NRC  Reg.  Guide  1.60  response  spectra.]  The  question  at  hand  is 
whether  these  response  measures  are  in  fact  independent  of  peak  accelera¬ 
tion.  To  examine  this  question,  the  140  records  are  ranked  according  to 
their  actual  peak  acceleration,  and  the  (conditional)  mean  and  standard 
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deviation  of  each  response  measure  are  calculated  for  subsets  of  records 

whose  peak  accelerations  lie  within  a  relatively  narrow  range  of  values. 

Formally,  suppose  that  Y  denotes  a  response  quantity  of  interest  and 

that  a  denotes  the  peak  acceleration.  Does  the  conditional  distribution 
P 

of  Y  depend  on  a  ?  The  conditional  response  statistics  can  be  estimated 
P 

from  a  subset  of  records  with  maximum  accelerations  close  to  a  .  Subsets 

P 

are  chosen  so  that  half  of  the  records  have  peak  accelerations  less  than  or 
equal  to  a^,  and  the  other  half  peak  accelerations  that  exceed  a^.  The  size 
of  each  subset  of  records  (n  =  20)  was  selected  on  the  basis  of  a  brief  sen¬ 
sitivity  study  in  which  the  number  of  records  in  each  subset  was  varied. 

The  analysis  procedure  is  illustrated  in  Figures  4  through  7.  The 
response  parameter  is  the  acceleration  amplification  factor  for  an  elastic 
one-degree  system  with  natural  period  T^  =  1.0  sec.  and  damping  ratio 
g  =  0.02.  Fig.  4  shows  the  amplification  factors  of  all  140  records  as  a 
function  of  their  peak  accelerations.  For  the  same  one-degree  system,  the 
sample  mean  of  the  amplification  factor  based  on  "moving"  subsets  comprising 
10,  20  and  30  records  is  shown  in  Figures  5,  6  and  7,  respectively. 

Throughout  the  report,  the  statistics  given  are  the  mean  m,  the  stan¬ 
dard  deviation  o,  and  the  coefficient  of  variation  V  =  a/m.  These  are  di¬ 
rectly  useful  in  evaluating  exceedance  probabilities  if  the  random  quan¬ 
tity  at  hand,  say  X,  is  assumed  to  be  Gaussian.  However,  the  lognormal 
distribution  is  frequently  an  attractive  alternative.  In  that  case,  the 
mean  and  the  standard  deviation  of  i.og  X  are  needed  to  evaluate  the  exceed¬ 
ance  probabilities.  If  Y  =  Hog  X,  the  mean  and  the  standard  deviation  of 

mx  and  the  coefficient  of  variation  are  related  as  follows 
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Figure  A.  Acceleration  amplification  factor  SA/ap  for  an  elastic  one-degree 
system  with  natural  period  =  1.0  sec.  and  damping  ratio  6  =  0.02 
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Figure  5.  Mean  acceleration  amplification  factor  SA/ap  for  Tn  -  1.0  sec.  and 
6  =  0.02.  Each  point  is  based  on  a  subset  of  10  records 


oy2  =  *og(Vx2  +  1) 

“y  =  iog  mx  -  \  ay2 

Hence,  if  the  statistics  of  X  are  given  (as  they  are  in  this  report)  the 
corresponding  statistics  of  iog  X  can  be  easily  evaluated. 

Results 

32.  This  section  of  the  report  presents  the  major  findings  of  the  statis¬ 

tical  study  in  the  form  of  a  series  of  figures  showing  the  statistics  of 
various  performance  parameters  (i.e.,  mean  m,  standard  deviation  o,  or 
coefficient  of  variation  V)  versus  peak  acceleration. 

a)  Figures  8-1  through  8-6  show  statistics  (m  and  m+o)  of  the  ac¬ 
celeration  amplification  factor  S./a  for  2%  damping  for  6  different 

A  p 

natural  periods  (T  =  0.1,  0.2,  0,5,  1.0,  2.0  and  5.0  sec.). 


b)  Figures  9-1  through  9-6  show  statistics  (m  and  m+o)  of  the  ac¬ 


celeration  amplification  factor  S./a  for  5%  damping  for  6  different 

A  p 

natural  periods  (T  =  0.1,  0.2,  0.5,  1.0,  2.0  and  5.0  sec.). 


c)  Figures  10-1  through  10-4  show  statistics  (m  and  nri-c)  of  the 

inelastic  acceleration  amplification  factor  S.(y)/a  for  a  one- 

A  p 

degree  system  with  natural  period  =  0.5  sec.  and  damping  5%,  for 
different  specified  values  of  the  ductility  factor  y  =  2,  3,  4  and 
5,  respectively. 

d)  Figures  11-1  through  11-4  show  the  mean  inelastic  acceleration 
response  rat^o  SA/SA(y)  for  4  different  values  of  ductility  factor 

y  (y  =  2,  3,  4  and  5).  The  different  figures  correspond  to  different 
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MEAN  &  MEAN  +  SIGMA  FOR  RESPONSE/PEAK  ACC  2  MEAN  &  MEAN  +  SIGMA  FOR  RESPONSE/PEAK  ACC 


gure  8-1.  Mean  and  mean  plus  standard  deviation  of  the  acceleration 
amplification  factor  S^/a^  for  2%  damping  and  Tr  =  0.1  sec. 


Figure  8-2.  Mean  and  mean  plus  standard  deviation  of  the  acceleration 
amplification  factor  S./a  for  2%  damping  and  T  =0.2  sec. 
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Figure  9-1.  Mean  and  mean  plus  standard  deviation  of  the  acceleration 
amplification  factor  SA/ap  for  5%  damping  and  Tn  =  0.1  sec. 
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Figure  9-2.  Mean  and  mean  plus  standard  deviation  of  the  acceleration 
amplification  factor  S  /a  for  5%  damping  and  T  =0.2  sec. 
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Figure  9-5.  Mean  and  mean  plus  standard  deviation  of  the  acceleration 
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Figure  9-6.  Mean  and  mean  plus  standard  deviation  of  the  acceleration 
amplification  factor  S  /a  for  5%  damping  and  Tr  =  5.0  sec. 
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igure  10-1.  Mean  and  mean  plus  standard  deviation  of  the  Inelastic 

cceleration  amplification  factor  S.(y)/a  for  5%  damping  and  natural 

A  p 

period  0.5  sec.,  and  for  a  ductility  factor  y  =  2 


Figure  10-2.  Mean  and  mean  plus  standard  deviation  of  the  inelastic 

acceleration  amplification  factor  S.(y)/a  for  5%  damping  and  natural 

A  p 

period  0.5  sec.,  and  for  a  ductility  factor  y  »  3 


1.1 


u 


PEAK  ACCELERATION  (G'S) 


Figure  10-3.  Mean  and  mean  plus  standard  deviation  of  the  inelastic 

acceleration  amplification  factor  S.(p)/a  for  5%  damping  and  natural 

A  p 

period  0.5  sec.,  and  for  a  ductility  factor  u  =  4 
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Figure  10-4.  Mean  and  mean  plus  standard  deviation  of  the  inelastic 

acceleration  amplification  factor  S. (p)/a  for  5%  damping  and  natural 

a  p 

period  0.5  sec.,  and  for  a  ductility  factor  u  =  5 


Figure  11-1.  Mean  inelastic  acceleration  response  ratio  sA/sA(u)  for 

four  different  values  of  the  ductility  factor  u  (p  =  2,  3,  4  and  5). 
The  damping  ratio  is  2%  and  the  natural  period  Tr  =  0.2  sec. 
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Figure  11-2.  Mean  inelastic  acceleration  response  ratio  sA/SA(p)  for 

four  different  values  of  the  ductility  factor  u  (u  =  2,  3,  4  and  5). 
The  damping  ratio  is  2%  and  the  natural  period  T  =0.5  sec. 
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Figure  11-3.  Mean  inelastic  acceleration  response  ratio  SA/SA(p)  for 

four  different  values  of  the  ductility  factor  u  (u  =  2,  3,  4  and  5). 
The  damping  ratio  is  2%  and  the  natural  period  Tn  =  1.0  sec. 
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PEAK  ACCELERATION  (G'S) 

Figure  11-4.  Mean  inelastic  acceleration  response  ratio  SA/SA(u)  for 

four  different  values  of  the  ductility  factor  u  (u  •  2,  3,  4  and  5). 
The  damping  ratio  is  2%  and  the  natural  period  Tfl  -  5.0  sec. 
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natural  periods  (T  =  0.2,  0.5,  1.0  and  2.0  sec.),  and  all  have  the 
same  2%  damping  ratio. 


e)  Figures  12-1  through  12-4  show  the  mean  inelastic  acceleration 

response  ration  S^/S^(y)  for  4  different  values  of  the  ductility 

factor  y  (y  =  2,  3,  4  and  5).  The  different  figures  correspond 

to  different  natural  periods  (T  =  0.2,  0.5,  1.0  and  2.0  sec.), 

n 

and  all  have  the  same  5%  damping  ratio. 


33.  The  next  series  of  figures  relate  to  the  variability  of  the  ac¬ 

celeration  amplification  factors,  S,/a  and  S.(y)/a  ,  for  simple  elastic 

A  p  A  p 

and  inelastic  one-degree  systems,  respectively: 

a)  Fig.  13  shows  the  coefficient  of  variation  V  of  the  accelera¬ 
tion  amplification  factor  S./a  for  2%  damping  for  3  different  na- 

A  p 

tural  periods  (T  =  0.2,  1.0  and  5.0  sec.). 

b)  Figures  14-1  through  14-4  show  the  standard  deviations  a  of  the 

elastic  amplification  factor  S./a  (top  curve  labeled  y=l)  and  the 

a  p 

inelastic  amplification  factor  S.(y)/a  (curves  labeled  y=2,  y=3, 

A  p 

etc.).  The  different  figures  correspond  to  different  natural  periods 

(T  =  0.2,  0.5,  1.0  and  5.0  sec.),  and  all  have  the  same  2%  damping 
n 

ratio. 


c)  Figures  15-1  through  15-4  show  the  standard  deviations  of  the 


elastic  amplification  factor  S,/a  (top  curve  labeled  y=l)  and 

A  p 


inelastic  amplification  factor  S.(y)/a  (curves  labeled  y=2,  v=3, 

A  p 


etc.).  The  different  figures  correspond  to  different  natural  periods 


(T  =  0.2,  0.5,  1.0  and  5.0  sec.),  and  all  have  the  same  5%  damping 


ratio. 
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Figure  12-1.  Mean  inelastic  acceleration  response  ratio  sA/sA(u)  for 

four  different  values  of  the  ductility  factor  p  (p  =  2,  3,  4  and  5). 
The  damping  ratio  is  5%  and  the  natural  period  is  Tr  =  0.2  sec. 
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PEAK  ACCELERATION  (G'S) 


Figure  12-2.  Mean  inelastic  acceleration  response  ratio  SA/SA(p)  for 

four  different  values  of  the  ductility  factor  p  (p  -  2,  3,  4  and  5). 
The  damping  ratio  is  5%  and  the  natural  period  T^  =  0.5  sec. 
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MEAN  FOR  INELASTIC  ACCELERATION  RATIO 


PEAK  ACCELERATION  (G'S) 

Figure  12-3.  Mean  inelastic  acceleration  response  ratio  SA/SA(p)  for 

four  different  values  of  the  ductility  factor  p  (p  =  2,  3,  4  and  5). 
The  damping  ratio  is  5%  and  the  natural  period  T^  -  1.0  sec. 
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Figure  12-4.  Mean  inelastic  acceleration  response  ratio  sa/sa(f)  for 

four  different  values  of  the  ductility  factor  p  (p  =  2,  3,  4  and  5). 
The  damping  ratio  is  5%  and  the  natural  period  T^  *  5.0  sec. 
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ure  15-1.  Standard  deviations  of  elastic  (y  =  1)  and  inelastic  (y  -  2, 
3,  4  and  5)  amplification  factors  for  5%  damping  and  =  0.2  sec. 


Figure  15-2.  Standard  deviations  of  elastic  (y  *  1)  and  inelastic  (y  =  2, 
3,  4  and  5)  amplification  factors  for  5%  damping  and  Tr  =  0.5  sec. 
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d)  Figures  16-1  and  16-2  show  the  coefficient  of  variation  V  of 
the  inelastic  acceleration  response  ratio  (i.e.,  ratio  of  elastic 
response  to  inelastic  response)  for  different  ductilities,  p=2, 
y=3,  etc.  The  two  figures  correspond  to  different  damping 

values  (6  =  0.02  and  6  =  0.05),  but  they  are  for  the  same  natural 
period  =  0.5  sec. 

The  final  series.  Figures  17-1  through  17-4, relate  to  the  assessment 
of  liquefaction  potential.  Each  figure  shows  a  scattergram  of  the  com¬ 
puted  "number  of  equivalent  cycles"  versus  the  actual  peak  ground  acceler¬ 
ation  for  140  accelerograms.  The  equivalence  is  based  on  different  pre¬ 
scribed  fractions  of  the  peak  acceleration  in  the  different  figures.  The 
fractions  for  which  results  are  provided:  0.55,  0.65,  0.75  and  0.85. 

Conclusions 

The  results  just  presented  provide  specific  quantitative  information 
about  the  effect  of  scaling  of  earthquake  records  on  seismic  performance 
predictions.  To  the  extent  that  the  set  of  records  used  in  the  analysis 
fairly  represents  the  body  of  information  about  measured  strong  earthquake 
ground  motion,  it  can  be  concluded  that  scaling  of  accelerogram  does  tend 
to  bias  performance  predictions.  If  normalized  w.r.t.  maximum  acceleration, 
the  response  of  simple  elastic  and  elasto-plastic  systems  and  the  equiva¬ 
lent  number  of  cycles  all  show  varying  degrees  of  dependence  on  the  maxi¬ 
mum  acceleration.  The  amount  of  error  introduced  by  scaling  varies  depend¬ 
ing  on  the  type  of  structure  (e.g.,  stiff  vs  flexible)  and  the  type  of 
seismic  performance  under  investigation.  For  all  but  very  stiff  struc¬ 
tures,  acceleration  amplification  factors  tend  to  decrease  as  the  actual 
peak  acceleration  increases. 
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Figure  16-1.  Coefficients  of  variation  of  the  inelastic  (p  ■ 
acceleration  amplification  factor  for  2%  damping  and  Tn 
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Figure  16-2.  Coefficients  of  variation  of  the  inelastic  (p  -  2,  3,  A  and  5) 
acceleration  amplification  factor  for  52  damping  and  T  =  0.5  sec. 
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number  of  cycles  versus  the  actual  peak  number  of  cycles  versus  the  actual  peal 

ground  acceleration  for  the  set  of  140  ac-  ground  acceleration  for  the  set  of  140 

celerograms.  The  equivalence  is  based  on  celerograms.  The  equivalence  Is  based 

the  acceleration  level  0.55  a  the  acceleration  level  0.65  a 
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The  results  confirm  the  validity  of  the  Newmark-Hall  procedure  for 
constructing  inelastic  response  spectra  by  multiplying  (elastic)  response 
spectra  by  an  inelastic  response  ratio.  The  inelastic  acceleration  response 
ratios  are  found  to  be  remarkably  independent  of  the  maximum  acceleration. 
However,  the  results  presented  here  do  provide  a  basis  for  suggesting  new, 
more  consistent  values  for  the  inelastic  acceleration  response  ratios. 

The  standard  deviation  or  coefficient  of  variation  (ratio  of  standard 
deviation  to  mean)  is  provided  for  each  seismic  response  parameter  studied. 
This  information  on  variability  is  useful  in  assessing  the  degree  of  con¬ 
servativeness  inherent  in (proposed) design  criteria  (whether  they  be  in 
terms  of  smooth  response  spectra  or  in  terms  of  selected  "design"  time 
histories) . 

It  is  not  suggested  that  the  dependence  of  the  various  (normalized) 
performance  measures  on  maximum  acceleration  is  direct  or  causative.  On 
the  contrary,  it  appears  that  much  of  the  variation  of  the  mean  (normalized) 
responses  in  function  of  maximum  acceleration  can  be  explained  in  terms 
of  the  strong-motion  duration  and  the  predominant  frequency.  These  ground 
motion  parameters  are  indirectly  related  to  peak  acceleration  through  their 
common  dependence  on  earthquake  magnitude,  epicentral  distance,  and  site 

ground  condition.  These  effects  are  further  explored  in  Part  IV  of  the 
report. 


PART  IV:  INTERPRETATION  OF  THE  RESULTS  OF 


THE  STUDY  ON  ACCELEROGRAM  SCALING 

Influence  of  Duration  and  Frequency  Content 

40.  A  second  statistical  study  was  carried  out,  based  on  the  140  accelero 
grams,  to  relate  the  strong-motion  duration  to  peak  acceleration  as  well 
as  to  magnitude  and  epicentral  distance.  Parameters  which  characterize 
the  frequency  content  of  the  ground  motion  were  analyzed  in  a  parallel  way. 

41.  In  earlier  reports  in  this  series  (1,  9,  16),  duration  was  taken  as 
the  "bracketed"  time  interval  in  which  acceleration  is  greater  than  0.05g. 
These  durations  are  listed  in  Table  1,  together  with  the  strong-motion 
durations  based  on  the  definition  proposed  by  Vanmarcke  and  Lai  (17) .  Ac¬ 
cording  to  this  definition,  duration  is  related  to  the  peak  acceleration 

a  =  max  I  a ( t) I  as  well  as  to  the  Arias  Intensity  I  =  fto  a2(t)dt.  The 
P  °  0 

resulting  duration  is  essentially  proportional  to  the  factor  . 

(Both  I  and  a  are  listed  in  Table  I).  Fig.  18  shows  a  scattergram  of  the 
o  p 

durations  according  to  the  two  definitions  for  those  accelerograms  for 
which  the  "bracketed  duration"  exists  (i.e.,  when  a^  0.05g) .  The  cor¬ 
relations  presented  below  are  based  on  the  Vanmarcke-Lai  definition  of  du¬ 
ration. 

42.  The  mean  duration  for  all  records  is  9.27  sec.,  and  the  standard 
deviation  is  8.7  sec.  For  the  32  near-field  records,  the  mean  and  the 
standard  deviation  are  6.34  and  5.55  sec.,  respectively.  For  the  far- 
field  records,  the  mean  is  10.2  sec.,  and  the  standard  deviation  is  9.4 
sec.  For  the  records  on  "rock,"  the  mean  duration  is  4.7  sec.,  and  for 
the  records  on  "soil,"  it  is  10.1  sec. 
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BRACKETED  DURATION 


Figure  18.  Scattergram  of  durations  based  on  two  different  definitions, 
the  "bracketed  duration"  (which  exists  for  records  with  ap  >0.05g)  and 
the  strong-motion  duration  proposed  by  Vanmarcke  and  Lai  (17) 
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43 


44. 


45. 


46. 


Although  there  is  considerable  scatter,  large  peak  accelerations 
tend  to  be  associated  with  small  strong-motion  durations.  Fig.  19  shows 
the  mean,  and  the  mean  plus  one  standard  deviation,  of  the  strong-motion 
duration  as  a  function  of  the  peak  acceleration.  As  before,  the  point 
on  the  "mean"  curve  corresponding  to  the  peak  acceleration  a^  is  obtained 
by  averaging  the  durations  for  a  subset  of  20  records  whose  "median"  peak 
acceleration  equals  a^.  Note  that  the  mean  duration  decreases  from  18 
sec.  to  about  3  sec.  as  the  (median)  value  of  ap  increases  from  0.015g 
to  0.25g.  The  coefficient  of  variation  (ratio  of  standard  deviation  to 
mean)  does  not  change  much  with  peak  acceleration  (see  Fig.  20). 

Scattergrams  of  the  "equivalent  number  of  cycles"  and  the  strong-mo¬ 
tion  duration  for  the  140  accelerograms  shown  in  Figures  21  through  24 
suggest  that  these  quantities  are  strongly  correlated.  These  scattergrams 
correspond  to  different  fractions  of  the  peak  acceleration,  ranging  from 
0.55  a  to  0.85  a  . 


P  P 

We  now  turn  to  the  dependence  of  duration  on  distance  and  magnitude. 
Fig.  25  shows  that  the  duration  tends  to  increase  with  distance  R.  The 
points  in  this  plot  again  correspond  to  the  statistics  (mean  or  mean  plus 
one  standard  deviation)  of  a  "moving  subset"  of  20  accelerograms  whose 
actual  distances  are  clustered  around  a  target  distance.  Fig.  26  indicates 
that  the  mean  duration  also  tends  to  increase  as  the  magnitude  M  increases, 
but  in  a  much  less  pronounced  way  (than  for  R) .  Figures  27  and  28  show 
the  dependence  of  duration  on  distance  for  records  from  earthquakes  with 
M  >_  6  and  M  <  6,  respectively. 

The  predominant  frequencies  of  the  set  of  140  accelerograms  do  not 
display  any  clear  dependence  on  peak  acceleration  (Fig.  29) .  There  is 
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EPICENTRAL  DISTANCE  (KM)  M>6  10  PTS 

Figure  27.  Mean  and  mean  plus  standard  deviation  of  strong-motion  duration 
versus  epicentral  distance  for  records  from  events  with  Richter  magnitude 
M  t  6.  (Calculations  based  on  a  "moving"  subset  of  10  records) 


EPICENTRAL  DISTANCE  (KM)  M<6  10  PTS 

Figure  28.  Mean  and  mean  plus  standard  deviation  of  strong-motion  duration 
versus  epicentral  distance  for  records  from  events  with  Richter  magnitude 
M  <  6.  (Calculations  based  on  a  "moving"  subset  of  10  records) 


PEAK  GROUND  ACCELERATION  (G)  20  PTS 

Figure  29.  Kean  and  mean  plus  standard  deviation  of  predominant  frequency 
versus  peak  acceleration.  (Calculations  are  based  on  a  "moving"  subset  of 

20  records) 


EPICENTRAL  DISTANCE  (KM)  20  PTS 

Figure  30.  Mean  and  mean  plus  standard  deviation  of  the  predominant 
frequency  versus  distance.  (Calculations  based  on  a  "moving"  subset 

of  20  records) 


also  no  clear  trend  in  function  of  earthquake  magnitude.  Since  high  fre¬ 
quency  components  of  ground  motion  decay  relatively  rapidly  with  distance, 
one  expects  the  predominant  frequency  to  decrease  with  distance.  This  ex¬ 
pected  trend  is  observed  in  Fig.  30. 

Conclusions 

47.  For  the  set  of  accelerograms  used  in  this  study,  large  peak  accelera¬ 
tions  tend  to  be  associated  with  small  durations.  The  duration  in  turn  in¬ 
fluences  seismic  response.  This  negative  correlation  between  peak  accelera¬ 
tion  and  duration  appears  largely  attributable  to  the  fact  that  both  peak 
acceleration  and  duration  depend  on  epicentral  distance.  The  reader  is 
cautioned  not  to  infer  from  this,  however,  that  ground  motions  having  large 
peak  accelerations  will  necessarily  have  small  duration.  For  example,  ground 
motions  from  very  large  magnitude  events  (which  are  hardly  represented  in 
the  body  of  instrumental  strong-motion  data)  might  have  both  large  accelera¬ 
tions  and  large  durations.  It  follows  that,  in  that  case,  the  larger  re¬ 
sponse  spectra  amplification  factors  (shown  in  the  statistical  study  to  be 
associated  with  small  maximum  ground  accelerations  except  for  very  stiff 
structures)  might  then  be  appropriate  for  use  in  constructing  site  response 
spectra  for  this  seismic  condition  (i.e.,  very  large  magnitude,  moderate 
distance) . 

48.  The  picture  that  emerges  is  that  duration  is  an  important  ground  motion 
parameter  which  is  presently  not  accounted  for  when  standard  response  spec¬ 
tra  shapes  are  used  in  the  specification  of  seismic  design  criteria  at  a  site. 
While  duration  has  a  pronounced  effect  on  response  spectra  amplification 
factors  (especially  for  moderate  and  long  periods,  or  when  damping  is  small), 
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it  will  often  have  a  decisive  effect  on  measures  of  response  of  systems 
susceptible  to  liquefaction,  low-cycle  fatigue,  accumulation  of  permanent 
deformation,  stiffness  degradation,  progressive  failure,  etc. 

The  observed  trends  of  the  response  amplification  factor  in  function 
of  peak  acceleration  can  also  be  explained  in  part  in  terms  of  the  indirect 
correlation  that  exists  between  the  predominant  frequency  of  a  ground  motion 
and  its  peak  acceleration.  For  an  earthquake  of  given  magnitude,  both  the 
peak  acceleration  and  the  predominant  frequency  tend  to  decrease  when  the 
epicentral  distance  increases.  This  makes  it  appear  that  the  predominant 
frequency  increases,  i.e. ,  that  there  is  relatively  more  high-frequency  con¬ 
tent  in  the  ground  motion,  when  the  peak  acceleration  increases.  This  ex¬ 
plains  why  the  spectral  amplification  factor  tends  to  increase  with  a^  for 
very  stiff  structures,  while  it  decreases  with  a^  for  long-period  structures. 

In  summary,  the  practice  of  using  standard  response  spectra  shapes 
scaled  on  peak  acceleration  fails  to  account  for  other  important  ground 
motion  parameters,  such  as  duration  and  predominant  frequency,  that  influ¬ 
ence  response  spectra  shapes.  By  the  same  token,  scaling  of  accelerograms 
and  response  spectra  on  peak  acceleration  introduces  systematic  errors  at¬ 
tributable  to  indirect  correlations  between  the  peak  acceleration  and  these 
other  ground  motion  parameters.  Estimates  of  these  errors  may  be  obtained 
from  the  plots  presented  in  Part  III.  The  sizes  of  these  errors  and  the  per- 
missable  values  of  scale  factors  will  vary  depending  on  the  seismic  response 
quantity  under  study  and  on  the  range  of  peak  accelerations  involved.  For 
example,  if  the  time  histories  are  used  for  liquefaction  studies,  scale 
factors  applied  to  accelerograms  should  perhaps  not  exceed  2,  while  for 
seismic  investigation  of  linear  elastic  systems,  factors  as  large  as  4  might 
be  tolerated. 
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PART  V:  COMPATIBILITY  BETWEEN  TIME  HISTORIES  AND 


SMOOTH  RESPONSE  SPECTRA 

51.  Compatibility  between  alternative  representations  of  earthquake  ground 
motion  at  a  site  is  an  important  practical  concern  in  earthquake  engineering. 
Consider  the  following  situation.  A  set  of  earthquake  records  has  been  se¬ 
lected  to  represent  a  prescribed  seismic  condition  at  a  site.  These  records 
will  be  used  as  input  for  time  history  analyses  of  seismic  response.  At 

the  same  site,  smooth  response  spectra  representing  essentially  the  same 
seismic  condition  are  also  needed  to  permit  response  spectrum-based  analy¬ 
sis.  The  question  is  how  one  should  proceed  to  derive  a  set  of  smooth  re¬ 
sponse  spectra  that  are  compatible  with  the  seismic  criteria  and  the  selected 
accelerograms.  It  is  assumed  that  the  time  histories  as  well  as  the 
response  spectra  will  be  scaled  with  respect  to  a  prescribed  peak  ground  ac¬ 
celeration. 

52.  It  is  common  in  present  practice  to  make  use  of  standard  response  spec¬ 
tra  shapes,  derived  on  the  basis  of  statistical  analysis  of  many  strong- 
motion  records  representing  a  wide  range  of  source  and  site  conditions  and 
epicentral  distances.  The  standard  response  spectra  shapes  may  correspond 
to  the  mean  of  the  amplification  factors  or  to  the  mean  plus  one  standard 
deviation.  Which  of  these  is  used  depends  largely  on  whether  or  not  it  is 
desirable  to  add  conservatism  at  this  particular  step  in  the  seismic  design 
process.  An  argument  can  also  be  made  for  always  using  mean  (or  median) 
response  spectra  and  scaling  the  spectra  upward  if  greater  conservatism  is 
sought . 

53.  There  is  no  unique  "best  way"  to  construct  smooth  response  spectra  that 
represent  a  set  of  time  histories  (scaled  to  peak  acceleration) .  It  is  as 
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much  an  art  as  it  is  a  science.  Obvious  factors  to  be  considered  are:  (a) 
the  size  and  quality  of  the  (recorded  ground  motion)  data  base,  (b)  desired 
simplicity  and  uniformity  of  the  response  spectra  shapes,  and  (c)  accuracy 
in  representing  expected  ground  motion  spectral  content.  Statistical  uncer¬ 
tainty  decreases  in  function  of  the  number  of  time  histories  that  are  used 
to  estimate  the  mean  and  the  standard  deviation  of  the  response  spectra  ordi¬ 
nates.  But  systematic  errors  or  biases  are  introduced  when  the  set  of  time 
histories  contains  records  with  source  or  site  characteristics  that  deviate 
significantly  from  the  seismic  condition  the  response  spectra  are  supposed 
to  represent. 

54.  Much  the  same  argument  can  be  made  when  one  analyzes  the  decision  about 
how  to  subdivide  the  frequency  range  for  the  purpose  of  computing  means  and 
standard  deviations  of  spectral  amplification  factors.  In  this  case,  the 
number  of  choices  can  be  limited  considerably  by  making  use  of  a  strong  body 
of  "prior  information"  that  indicates  that  response  spectra  are  approximate¬ 
ly  proportional  to  ground  acceleration,  ground  velocity  and  ground  displace¬ 
ment  in  different,  adjacent  frequency  ranges.  In  other  words,  it  may  be  as¬ 
sumed  at  the  outset  that  the  smooth  response  spectra  will  be  piecewise  linear , 
and  alternately  parallel  to  the  constant  acceleration,  velocity  and  displace¬ 
ment  lines  in  the  usual  tripartite  logarithmic  representation  of  response 
spectra.  It  will  therefore  suffice  to  calculate  (for  each  damping)  the  mean 
and  the  standard  deviation  of  a  single  amplification  factor  that  can  be  ap¬ 
plied  throughout  the  frequency  range  where  (that  factor  is)  essentially  con¬ 
stant. 

55.  The  proposed  methodology  is  introduced  through  an  example  involving 
two  hypothetical  seismic  design  situations.  Each  site/source  condition  is 
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specified  in  terms  of  magnitude,  epicentral  distance,  and  local  site  con¬ 
dition  (soil  vs  rock) .  Also,  based  on  information  about  motion  attenua¬ 
tion  with  distance  and  dependence  on  magnitude,  a  "design"  maximum  ground 
acceleration  is  specified  for  each  case.  The  criteria  for  the  (hypotheti¬ 
cal)  seismic  conditions  are  summarized  in  Table  2.  In  the  first  case,  a 
moderate  magnitude  -  short  distance  event  is  specified,  while  the  second 
case  involves  a  large  magnitude  -  long  distance  event.  The  particular  at¬ 
tenuation  relations  used  in  Table  2  are  based  on  western  United  States 
data.  It  must  be  emphasized  that  they  serve  only  an  illustrative  purpose 
here.  Detailed  discussion  about  methodology  to  obtain  ground  motion  ampli¬ 
tudes  given  a  seismic  condition,  or  to  obtain  design  amplitudes,  is  beyond 
the  scope  of  this  report.  A  pertinent  reference  is  the  report  by  Krinitz- 
sky  and  Chang  in  this  series  (1) .  It  summarizes  information  about  the  re¬ 
lations  between  epicentral  M.M.  Intensity  and  magnitude,  between  M.M.  In¬ 
tensity  and  distance,  and  between  maximum  ground  motion  amplitudes  and 
M.M.  Intensity.  The  report  emphasizes  the  distinction  that  exists  between 
the  near  field  and  the  far  field  as  regards  the  relation  between  motion 
amplitudes  and  M.M.  Intensity.  It  presents  upper  bounds  of  observed  ground 
motion  amplitudes  as  well  as  estimates  of  the  conditional  probability  of 
exceeding  a  motion  amplitude  given  the  M.M.  Intensity.  Needless  to  say, 
this  is  precisely  the  kind  of  information  the  engineer  needs  to  make  an 
assessment  of  site  ground  motion  amplitudes  during  a  specified  seismic 
condition. 
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Table  2 


Criteria  for  Two  Hypothetical  Seismic 
Design  Conditions 


Site/Source  Conditions 

Case  1 

Case  II 

Magnitude 

6.5 

8.0 

Focal  distance,  km 

50 

150 

Local  site  geology 

Soil 

Soil 

Maximum  acceleration* 

(1) 

99.5 

69.9 

predicted,  cm/sec^ 

(2) 

118 

102 

For  design 

O.lg 

O.lg 

*  Maximum  accelerations  are  predicted  from  attenuation  rela¬ 
tionships  proposed  by  Donovan  (1973)  and  Orphal  and  Lahoud 
(1974).  These  relationships  all  have  the  form: 

,  boM  D  -bo 
a  =  b,  e  ‘  R  J 
pi  o 

2 

where  a  =  maximum  acceleration  in  cm/sec  ,  M  =  Richter  Magni- 
P 

tude,  Rq  =  R  +  c,  R  =  focal  distance  (km),  b^,  b^,  b^  and  c  are 

constants.  The  respective  values  of  the  constants  (b^,  b^,  b^ 

and  c)  are:  (1)  Donovan:  1080,  0.51,  1.32  and  25;  and  (2) 
Orphal  and  Lahoud:  69,  0.92,  1.39  and  0. 
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It  is  assumed  that  the  pool  of  available  earthquake  records  consists 
of  140  horizontal  components  of  strong-motion  records,  the  same  data  set 
used  in  the  study  of  the  effect  of  accelerogram  scaling.  The  records  se¬ 
lected  to  represent  the  site  ground  motion  are  given  in  Tables  3  and  4. 
Obviously,  the  selection  process  requires  compromise  and  judgement,  but 
the  following  objective  criteria  were  used:  (a)  agreement  of  site  geology 
(rock  vs  soil) ,  (b)  close  matching  of  the  specified  magnitude  and  distance 
and  (c)  minimizing  the  amount  of  accelerogram  scaling  required.  In  gene¬ 
ral,  other  constraints  may  have  to  be  added,  e.g.,  the  type  of  faulting, 
focal  depth,  regional  geology.  Of  course,  the  more  constrained  the  search 
the  harder  it  will  be  to  find  the  desired  number  of  motions  without  devi¬ 
ating  too  much  from  the  specified  site  condition.  In  the  cases  at  hand,  8 
accelerograms  (2x4  horizontal  components)  were  selected  to  represent  the 
moderate  magnitude  -  short  distance  event  (Case  I).  For  the  large  magni¬ 
tude  -  large  distance  event  (Case  II)  only  4  accelerograms  (2x2  horizon¬ 
tal  components)  were  found  that  came  reasonably  close  to  the  prescribed 
site/source  condition.  Tables  3  and  4  list,  for  each  accelerogram  selec¬ 
ted,  the  peak  acceleration  and  velocity  (a^  and  v^) ,  event  magnitude, 
source- to-site  distance,  the  site  condition,  and  the  scale  factor  required 


Table 


Average:  0.11  g  Mean:  129.8 

St.  Dev:  60.11 
C.o.v.:  0.463 


the  "Case  II"  Site  Seismic  Condition 


to  match  the  "design"  maximum  acceleration. 

57.  The  next  step  is  to  assess  the  implications  of  the  selection  of  ground 

motions  on  the  response  spectra  shapes.  The  response  spectra  are  evaluated 

for  2%  and  5%  damping  and  for  the  12  natural  periods  ranging  between  0.05 

and  2  sec.  The  corresponding  spectral  amplification  factors  S ,/a  and  S  /v 

A  p  v  p 

are  listed  in  Tables  5  through  8.  The  velocity  amplification  factor  S^/v ^ 
is  directly  related  to  the  acceleration  amplification  factor  S^/a^  as  follows: 


v  2tt  v  a 
P  P  P 


58.  The  mean  amplification  factors  are  listed  just  below  each  column  of 
amplification  factors  in  Tables  5  through  8.  The  corresponding  standard 
deviations  and  coefficients  of  variation  are  listed  below  the  mean  values. 
It  will  now  be  shown  how  the  data  presented  in  these  tables  can  be  used  to 
construct  smooth  response  spectra.  Consider  as  a  specific  example  the  data 
in  Table  5  which  pertain  to  the  2%  damped  response  spectra  for  Case  I. 

59.  In  Table  5,  consider  first  the  values  of  the  coefficient  of  variation 

(c.o.v. )  of  the  ratios  S./a  and  S,,/v  as  a  function  of  natural  period. 

A  p  V  p 

The  c.o.v.  values  of  the  velocity  amplification  factor  are  significantly 
lower  within  a  range  of  moderate  periods  (from  about  0.5  to  1.0  sec.)  than 
on  either  side  of  this  period  range.  The  c.o.v.  values  of  the  acceleration 
amplification  factor  do  not  vary  much  between  0.1  and  0.4  sec.  but  make  an 
upward  jump  (at  T^  =  0.5  sec)  just  when  the  variability  of  the  velocity  am¬ 
plification  factor  begins  to  decrease  dramatically.  Based  on  these  patterns 
(of  dependence  of  c.o.v.  on  natural  period),  the  "corner  period,"  which  will 
separate  the  period  ranges  of  amplified  acceleration  and  amplified  velocity. 
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can  be  roughly  estimated  (the  corner  period  is  about  0.5  sec.  in  this 
case) . 

60.  The  second  step  is  to  evaluate  the  mean  and  the  standard  deviation 
of  the  amplification  factors  S^/a^  and  S^/v^.  Note  the  patterns  of  vari¬ 
ation  of  the  mean  amplification  factors  with  period.  The  mean  values  of 

S./a  initially  increase  with  period,  then  reach  a  plateau  in  the  range 
A  p 

0.2  sec.  <_  Tr  <_ 0.5  sec.,  and  finally  decrease  rapidly  as  increases  be¬ 
yond  0.5  sec.  In  the  range  0.2  sec.  <  T  <  0.5  sec.,  the  mean  of  the  ac- 

—  n  — 

celeration  amplification  factor  is  2.9  and  its  standard  deviation  is  about 
0.8.  Similarly,  in  the  amplified  velocity  range,  0.5  sec.  <_  T^  <_  2  sec., 
the  mean  of  the  velocity  amplification  factor  is  2.05  and  its  standard 
deviation  is  about  0.8. 

61.  In  the  short  period  range  of  the  response  spectrum,  a  gradual  transi¬ 
tion  must  be  made  between  the  zero-period  value  S . /a  =1  and  the  level 

A  p 

prevailing  in  the  period  range  where  the  acceleration  amplification  is  con¬ 
stant.  A  solid  guide  is  the  array  of  computed  mean  acceleration  amplifi¬ 
cation  factors  in  this  transition  range  (see  Table  5) .  The  corresponding 
standard  deviations  are  also  available  to  aid  in  constructing  scaled  spec¬ 
tra  corresponding  to  the  mean  plus  one  standard  deviation. 

62.  The  following  procedure  is  suggested  for  estimating  the  mean  accelera¬ 
tion  amplification  factor  in  the  amplified  velocity  range  of  the  response 

spectrum:  (a)  estimate  the  ratio  v  / a  either  by  averaging  the  values  of 

P  P 

this  ratio  for  each  of  the  records  selected  or  by  using  the  attenuation 

relationships  (predictions  of  a  and  v  in  function  of  M  and  R) ,  and  (b) 

P  P 

multiply  this  ratio  by  the  mean  amplification  factor  S^/v^,  and  (c) 

multiply  the  result  by  w  =  (2tt/T  ) .  A  parallel  procedure  may  be  followed 

n  n 

in  the  amplified  displacement  range  of  the  response  spectrum. 
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63.  It  is  proposed  that  the  uncertainty  in  the  ratio  vp/ap  n°t  be  re¬ 
flected  in  the  standard  deviation  of  the  scaled  response  spectra  in  the 
amplified  velocity  range.  The  standard  deviation  need  only  reflect  uncertain¬ 
ty  in  the  amplification  factor  S„/v  .  This  stems  from  the  fact  that  the 

V  p 

maximum  ground  velocity  can  be  estimated  independently,  just  as  well  as 
the  maximum  ground  acceleration.  In  the  velocity  range  of  the  response 
spectrum,  when  it  comes  to  assessing  variability,  one  may  view  the  ground 
velocity,  not  the  ground  acceleration,  as  the  "anchor."  The  result  is  that 
the  coefficient  of  variation  of  the  scaled  response  spectra  will  not  nec¬ 
essarily  be  larger  for  moderate  and  large  periods  than  for  relatively  short 
periods. 

64.  The  smooth  response  spectra  for  2%  and  5%  damping  (for  natural  periods 
below  2  sec.)  corresponding  to  the  two  seismic  conditions  are  shown  in 
Figures  31  and  32,  respectively.  Each  figure  shows  the  response  spectra 
corresponding  to  the  mean  and  to  the  mean  plus  one  standard  deviation.  A 
direct  comparison  of  the  mean  5%  damped  response  spectra  corresponding  to 
the  two  seismic  conditions  is  made  in  Fig.  33. 

65  We  have  just  outlined  how  "compatible"  smooth  response  spectra  can 

be  estimated  from  a  set  of  time  histories  representing  a  well-oefined 
seismic  condition.  For  the  sake  of  completeness,  it  should  be  mentioned 
that  it  is  also  possible  to  obtain  compatible  seismic  input  in  the  form 
of  the  spectral  density  function  and  strong-motion  duration  ,  which  con¬ 
stitute  the  input  for  seismic  random  vibration  analysis  of  structures. 

The  duration  and  the  parameters  of  the  spectral  density  function  may  now 
be  viewed  as  random  variables  whose  properties  (mean,  standard  deviation) 
can  be  estimated  partly  from  correlations  with  magnitude  and  distance,  and 
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PSEUDO-VELOCITY  RESPONSE  SPECTRUM,  IN./SEC. 


NATURAL  PERIOD,  SEC  . 

Figure  31.  Smooth  response  spectra  corresponding  to  mean  and  mean  plus 
standard  deviation  for  2%  and  5%  damping  for  "Case  I"  site  seismic 

condition 
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partly  from  the  specific  earthquake  records  selected  to  represent  the 
seismic  condition  under  study.  This  probabilistic  input  representation 
is  also  a  starting  point  for  generating  compatible  artificial  earthquake 
time  histories.  The  details  of  these  procedures  are  beyond  the  scope  of 
this  report. 

Design  Response  Spectra 

Throughout  Part  V  of  the  report,  reference  has  been  made  to  compatible 
ground  motion  representations  for  a  given  site  seismic  condition.  For  pur¬ 
poses  of  illustration,  two  different  seismic  conditions  (Cases  I  and  II) 
were  considered.  The  terminology  "design  seismic  condition"  or  "design  re¬ 
sponse  spectra"  was  deliberately  avoided.  The  smooth  response  spectra  cor¬ 
responding  to  a  given  seismic  condition  may  indeed  be  interpreted  as  design 
response  spectra  if  that  particular  condition  greatly  dominates  the  seismic 
risk  at  the  site.  If  two  or  more  different  seismic  conditions  significantly 
contribute  to  the  seismic  risk  at  the  site,  then  it  becomes  necessary,  in 
order  to  generate  proper  design  response  spectra,  to  combine  the  information 
about  the  (conditional)  responses  given  each  seismic  condition  with  informa¬ 
tion  (provided  by  seismic  risk  assessment)  about  the  relative  likelihood  of 
occurrence  of  each  seismic  condition.  A  procedure  to  do  this  is  illustrated 
below. 

Take  as  an  example  a  "soil  site"  where  the  0.1-g  design  earthquake  may 
originate  either  from  a  moderate  earthquake  at  moderate  distance  (Type  I)  or 
from  a  large,  distant  earthquake  (Type  II).  The  probabilities,  Pj  and  pIT  = 
1  -  Pj,  that  the  design  earthquake  will  be  of  a  particular  type  can  be  ob¬ 
tained  from  a  combination  of  seismological  and  geological  information,  sup- 
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ported  by  judgement  and  perhaps  by  formal  seismic  risk  analysis.  Ground 
motions  and  smooth  site  response  spectra  representing  the  two  seismic  con¬ 
ditions  have  been  obtained  before.  Specifically,  means  and  standard  devi¬ 
ations  of  the  amplification  factors  governing  the  shape  of  the  response 
spectra  in  the  various  period  ranges  have  been  determined.  From  this  in- 
fo mat ion ,  one  can  evaluate  the  probability  that  the  amplification  factor 
SA/an  is  less  than  a  specified  value,  say  r,  given  the  type  of  earthquake. 
These  conditional  probabilities  will  be  denoted  by  p(r|l)  and  p(rjll). 

The  probability  that  the  design  response  spectra  amplification  fac¬ 
tor  will  be  less  than  a  given  value  r  is 

P(r)  =  p(r|l)pI  +  pCrllDpjj  . 

This  probability  of  no  exceedance  must  be  evaluated  for  a  range  of  values 
of  the  factor  r.  From  the  relationship  between  p(r)  and  r,  the  engineer  can 
then  obtain  the  design  amplification  factor  corresponding  to  a  prescribed 
exceedance  probability.  To  carry  out  the  computations,  a  choice  has  to  be 
made  concerning  the  probability  distribution  of  the  amplification  factors 
for  each  type  of  earthquake.  This  choice  is  not  critical  if  probabilities 
of  no  exceedance  between  0.5  and  0.9  are  involved.  ^ The  normal  distribution 
was  assumed  in  the  illustrative  calculation.^  Of  course,  the  procedure  needs 
to  be  repeated  for  a  limited  number  of  periods  for  each  damping  ratio  for 
which  response  spectra  are  sought.  As  an  illustration,  the  median  design  re¬ 
sponse  spectra  were  calculated  for  a  situation  where  the  large,  distant  earth¬ 
quake  (Type  II)  was  judged  to  carry  30%,  and  the  closer,  moderate  earthquake 
(Type  I)  70%,  of  the  seismic  risk  at  the  0.1-g  design  level.  It  follows  that 
Pj  »  0. 7  and  p^  =  0.3.  The  resulting  median  design  response  spectra  for  5% 
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damping  are  shown  in  Fig.  33,  together  with  the  median  (or  mean)  response 
spectra  for  the  two  "constituent"  seismic  conditions. 

Conclusion 

69.  In  Part  V  of  the  report,  a  procedure  has  been  presented  for  developing 
site-specific  smooth  response  spectra  that  are  in  agreement  with  (and  are 
derived  from)  a  set  of  accelerograms  selected  on  the  basis  of  their  compati¬ 
bility  with  a  prescribed  site  seismic  condition.  Such  a  condition  may  be 
described  in  terms  of  an  event  magnitude  and  a  distance  (or  a  range  of  dis¬ 
tances)  as  well  as  in  terms  of  the  soil  conditions  at  the  site;  of  course, 
further  constraints  could  be  imposed,  e.g.,  on  the  source  mechanism  or  the 
tectonic  province.  An  alternative  is  to  express  the  site  seismic  condition 
in  terms  of  maximum  ground  motion  amplitudes  and  duration.  Clearly,  the  pro¬ 
cedure  permits  all  potentially  significant  factors  to  be  considered  in  con¬ 
structing  te-specific  spectra.  It  differs  in  this  sense  from  the  approach 
presented  by  Seed  et  al.  (18) ,  which  singles  out  local  site  geology  as  the 
factor  affecting  the  shape  of  site-specific  response  spectra. 

70.  The  proposed  approach  to  response  spectra  smoothing,  while  obviously  not 
unique,  has  the  following  desirable  features:  (a)  it  relies  on  "prior"  infor¬ 
mation  about  response  spectra  shapes  in  the  acceleration,  velocity  (and  dis¬ 
placement)  ranges  of  the  spectrum,  (b)  it  permits  the  calculation  of  rela¬ 
tively  stable  estimates  of  the  means  and  standard  deviations  based  on  only 

a  small  number  of  records,  and  (c)  it  makes  full  use  of  information  about 
maximum  ground  velocity  (or  displacement)  and  about  the  amplification  fac¬ 
tors  in  the  velocity  (or  displacement)  range  of  the  response  spectra.  Fi¬ 
nally,  the  report  outlines  the  steps  required  to  obtain  smooth  design  re¬ 
sponse  spectra  if  two  or  more  different  types  of  earthquakes  can  cause  de¬ 
sign  level  shaking  at  the  site. 
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